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Ground State Phase Diagram of 2D Electrons in High Magnetic Field 
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The ground state of 2D electrons in high magnetic field is studied by the density matrix 
renormalization group method. The ground state energy, excitation gap, and pair correlation 
functions are systematically calculated at various fillings in the lowest and the second lowest 
Landau levels. The ground state phase diagram, which consists of incompressible liquid state, 
compressible liquid state, stripe state, pairing state, and Wigner crystal is determined. 
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§1. Introduction 

Electrons in high magnetic field in two dimensional 
systems exhibit various interesting phenomena as a typ- 
ical many body interacting system. In two dimen- 
sional systems, kinetic energy of electrons is completely 
quenched by perpendicular magnetic field and remaining 
macroscopic degeneracy is lifted by Coulomb interaction. 
This is the origin of many exotic behaviors of quantum 
Hall systems and the reason why much effort has been 
needed to understand these systems. 

As typical ground states in strong magnetic field, 
Laughlin state and Wigner crystal have attracted much 
attention. Laughlin state is a quantum liquid with an 
excitation gap, which is an exact ground state for stroiig 
short range repulsive interaction at = l/(2m + 1)J3) 
On the other hand, Wigner crystal is the ground state 
of classical point particles realized in the limit of strong 
magnetic field. Their relative stability depends on the ra- 
tio between the magnetic length I = (Tic/ eBY^'^ ^ which 
is the length scale of one particle wave function, and the 
mean distance between the electrons ~ n~^/^ given by 
the density n of two-dimensional electrons. Thus the 
determination of the ground state phase diagram with 
respect to the filling factor v = 27r^^n is an interesting 
issue for understanding the competition between com- 
pletely different electronic states. 

Experimentally, clear fracticmal quantum Hall effect 
has been observed at u =J./3,B' which is naturally ex- 
plained by Laughlin state,EP and exact numerical diago- 
nalizations confirmed its existence in the lowest Landau 
level. The formation of Wigner crystal has been sup- 
ported by experiments on transport properties at low 
fillings u ^ 1/5, where observed insulating behavior is 
explained by the pinning of Wigner crystal by residual 
potential fluctuations.Q'B^ 

With decreasing magnetic field, part of electrons oc- 
cupy higher Landau levels. In higher Landau levels, 
one particle wave function extends over space with os- 
cillations, and this change modifies effective interaction 
between the electrons in each Landau level. Since ex- 



tended one particle wave function generates long range 
exchange interaction, Laughlin state stabilized by strong 
short range interaction will be unstable. Instead, long 
range interaction favors charge density waves (CDW's). 
Within a Hartree-Fock approximation, Koulakov et al. 
showed that CDW's called 'stripe' and 'bubble' are sta- 
ble against Laughlin state in higher Landau levels with 
index N > 2.13'LP After their calculations, anisotropic re- 
sistivity consistent with uniaxial electroniiistate of stripe 
order has been experimentally observed.Botj^ 

In the second lowest Landau level, more interesting 
ground states are expected. At vn = 1/2, where vn is 
the filling factor of partially filledJLpidau level, BCS type 
pairing state has been proposedtiHiS'tj^ to epqplain even 
denominator fractional quantum Hall effect .tS Around 
UN = 0.4 and — 0.24, ree^itrant integer quantum Hall 
effects have been observed.ll3) At un = 1/3 and-JZ5i, 
fractional quantum Hall effects have been observedJlil't3 
The origin of these phenomena is still unclear and sys- 
tematic study is needed. 

In order to identify diverse ground states and deter- 
mine the ground state phase diagram, we have to solve 
quantum many body problem of Coulomb interaction. In 
the present study we use the density matrix renormal- 
ization group (DMRG) algorithnitSllZP to deal with large 
size systems extending limitation of exact diagonaliza- 
tions, and systematically calculate the ground state wave 
function and excitation energies for the lowest and the 
second lowest Landau levels. Based on the pair correla- 
tion functions and the excitation energies, we determine 
the ground state phase diagram, which consists of incom- 
pressible liquid state, compressible liquid state, stripe 
state, paring state, and Wigner crystal. 

§2. Model and Method 

The Hamiltonian for electrons in Landau levels con- 
tains only Coulomb interaction. After the projection 
onto the Landau level with index the Coulomb in- 
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teraction is written as 

i<3 q 

(2.1) 

where is the guiding center coordinate of the 
iih electron, which satisfies the commutation relation, 
= il^Sjk^ Ln{x) are the Laguerre polynomi- 
als, and V{q) = 2i\e^ jeq is the Fourier transform of the 
Coulomb interaction. The magnetic length i is set to be 
1, and we take jet as units of energy scale. We omit 
the component at g = 0, which is canceled by uniform 
positive background charge. We assume completely spin 
polarized ground state and neglect the electrons in fully 
occupied Landau levels. We also assume that the width 
of the wave function perpendicular to the two dimen- 
sional plane is sufficiently small compared with i. 

In order to obtain the ground stale .of the Hamilto- 
nian, we employ the DMRG methodJi^'tlP This method 
enables us to calculate ground state wave function of 
large systems with controlled accuracy. The algorithm 
of this method is summarized as follows: We start from 
a small system, i.e. a system consisting of only four lo- 
cal orbit als. We divide the system into two blocks, and 
add new orbitals at the end of two blocks to expand the 
blocks. We then calculate the ground state wave func- 
tion of the total system and obtain the density matrix of 
each expanded blocks. We then restrict the basis states 
in the expanded blocks by keeping only eigenstates of 
large eigenvalues of the density matrix. We add new or- 
bitals again and repeat the above procedure until we get 
desired size of system. We then use the finite system al- 
gorithm of the DMRG to improve the ground state wave 
function. After several sweeps, we obtain the most op- 
timal ground state within the restricted number of basis 
states. The error due to the truncation of basis states 
is estimated from the eigenvalues of the density matrix, 
and the accuracy of the wave function is systematically 
improved by increasing the number of basis states kept 
in the blocks. 

§3. Results 

In the following, we calculate ground state wave func- 
tion of the Hamiltonian for Landau levels of = and 1. 
We study various size of systems with up to 25 electrons 
in the unit cell of x Ly with periodic boundary con- 
ditions in both X and y directions. We choose the aspect 
ratio Lx/Ly by the position of the energy minimum with 
respect to L^/Ly in order to avoid artificial determina- 
tion of the periodicity of CDW's. The truncation error in 
the DMRG calculation is typically 10~^ for 25 electrons 
with 180 states in each blocks. The existing results of 
exact diagonalizations are completely reproduced within 
the truncation error. Since the present Hamiltonian has 
the particle-hole symmetry, we only consider the case of 
vn < 1/2. 

3.1 Compressible liquid at v = 1/2 

We first investigate the ground state at = 1/2 in the 
lowest Landau level. Since the denominator of v is even 
and the filling factor is large, the ground state is neither 



Laughlin state nor Wigner crystal. Based on the compos- 
ite fermion theory, the flux attached to the electrons com- 
pletely cancels external magnetic field within the pj^an 
field level and Fermi liquid ground state is expected. Il3'tj) 
Indeed, experiments show satmration in longitudinal re- 
sistivity at low temper atures,E2P which is consistent with 
Fermi liquid like compressible ground state. 

In order to confirm the compressible ground state, 
we first analyze size dependence of the excitation gap. 
The obtained results for systems with up to 25 elec- 
trons are shown in Fig. 1. Although the excitation 
gap does not follow single scaling function, all the re- 
sults are smaller than the upper bound shown by the 

— 1/2 

dashed line, which is proportional to A's , where Ng 
is the number of one-particle states in the unit cell. 
We thus expect vanishing of the excitation gap in the 
thermodynamic limit. We find small excitation gap at 
small number of electrons. Such vanishing excitation 
gap appears for free electrons in the square unit cell at 
A^e = 2,3,4,6,7,8,10,11,12,14,15,..., with periodic bound- 
ary conditions. This means large excitation gap appears 
only at A'e = 1,5,9,13,21,25,..., where electrons occupy all 
the degenerate states forming closed shell. The present 
result of interacting system also shows a shell structure 
but it is not the same to that of free electrons in con- 
trast to the results in the spherical geometry, where the 
same shell structure is obtained.EiP Nevertheless, the es- 
timated effective mass 1/m* of the composite fermion 
determined from the value of the gap at A'g = 9 and 25 
is 0.12 and 0.27, respectively, in units of e^/e^, which are 
roughly consistent with the result 0.20 obtained from the 
exact, diagonalizations in the system of spherical geome- 
try.!^ 

In Fig. 2 (a) we show the ground state pair correlation 
function defined in the following equation: 

^^''^ ^ NAN J- 1) ^ + ~ (^-^^ 

where |^) is the ground state. The correlation function 
is circularly symmetric and short range power low expo- 
nent is two as shown in Fig. 3. These results are also 
consistent with gapless liquid ground state. 

3.2 Incompressible liquids at u = n/{l -\- 2n) 

When we decrease filling factor from 1/2, fractional 
quantum Hall effects ap. experimentally observed at var- 
ious fractional fillings .E3'E3) In particular, the primary 
series oi v = n/(l + 2n) with n = 1,2,3,4,5 is well 
known for its large excitation gap. Here, we confirm the 
existence of excitation gap at = n/(l + 2n) by system- 
atically calculating the lowest excitation gap for filling 
factors between u = 0.3 and 0.5. The obtained result for 
systems with 12, 14, 15, 16, 18, 20, 24, and 25 electrons 
are shown in Fig. 4. As is clearly seen, large excita- 
tion gap is obtained for u = n/(l + 2n). The size of 
the gap decreases with increasing n, and its dependence 
is consistent with experiments. Based on the composite 
fermion theory, the effective mass is estimated from the 
excitation gap. The estimated value of 1/m* in units 
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Fig. 1. Size dependence of the lowest excitation gap at u = 1/2 
for various number of electrons Ne in square unit cell. Ng is the 
number of one-particle states in the unit cell. The dashed line is 
a guide for eyes. 



of e'^ /s£ is around 0.23, which is almost the same to the 
value ~ 0.2 determined at u = 1/2. We mention that the 
composite fermion picture also prpiicts small excitation 
gap in finite system at = 4/ll,ElP which is consistent 
with the present result shown in Fig. 4. 

In order to confirm liquid ground state at v = n/{l -\- 
2n), we next see the ground state pair correlation func- 
tions. The results for u = 5/11,2/5, and 1/3 are shown in 
Figs. 2 (b), (c) and (f), respectively. As is seen, the corre- 
lation functions are circularly symmetric, which confirms 
liquid ground state. 
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Fig. 2. Ground state pair correlation functions in the lowest Lan- 
dau level for various u. Ne = 20 for u = 1/2, 5/11, 2/5, 10/27, 
5/14, and 1/3. A^e = 16 for u = 4/13 and 1/4. 



3.3 Weak stripes below v ~ 0.42 

We next consider almost gapless ground state between 
the incompressible liquid states. In Figs. 2 (d), (e), (g), 
and (h), the pair correlation functions at z^ = 10/27, 
5/14, 4/13 and 1/4 are shown. We find weak stripe cor- 
relations. Similar stripes are observed for v < 10/24 
except around the incompressible states at v = 2/b and 
1/3. The stripe correlations are gradually enhanced with 
decreasing v as shown in Fig. 5. We think the gradual 
enhancement in the stripe correlations are due to the 
instability to Wigner crystal, and quantum fiuctuations 
partially melt Wigner crystal to form stripe ground state. 
In fact, short range correlations below r ~ 5f are almost 
the same to that of Wigner crystal. This is clearly dif- 
ferent from the stripe state observed in higher Landau 
levels, where anisotropic resistivity is experimentally ob- 
served. The difference from the stripe state in higher 
Landau levels is clearly shown in §3.4. We note that 
the present results are different from those obtained in 
recent studies around v where stripe state with 

much longer period of 10^3^ and bubble state around 
y = 3/8E3' are predicted based on the composite fermion 
theory. The present DMRG results show that the period 




Fig. 3. Pair correlation functions in the lowest Landau level for 
various u near half filling. Inset shows logarithmic plot near the 
origin. 
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Fig. 4. The lowest excitation gap at various u. The number of 
electrons A^e is between 12 and 25. 



Fig. 6. Ground state pair correlation functions at u = 1/5 for 
various 5V^. Ne = 12. dV^, = corresponds to pure Coulomb 
interaction. 




Fig. 5. Ground state pair correlation functions in the lowest Lan- 
dau level at various h>. Ne = 24 for = 3/8. Ne = 16 for i/ = 
4/13 and 1/4. A^e = 12 for ly = 1/5 and 1/6. 



of the stripes at = 3/8 is about as shown in Fig. 3 
and the stripe state is stable in wide range below u ~ 0.42 
except around incompressible state at z/ = 2/5 and 1/3. 
Detailed study at z^ = 3/8 with small size dependence of 
the stripes will be published elsewhere. 

3.4 Laughlin state and stripe order at v = 1/5 

The stripe correlation in the ground state is enhanced 
with decreasing v. The correlation function is anisotropic 
even at z^ = 1/5 as shown in Fig. 6 (b), and it is not clear 
whether there exist transition to incompressible state 
around z/ ^ 1/5. In order to identify the ground state at 
V = 1/5, we introduceJHaldane's pseudo potentials Vm 
between the electrons.ElP It is shown that the Laughlin 
state at z^ = 1/5 is an exact ground state for strong short 
range repulsive interaction consisting of only Vi and V3, 
which act only between electron pairs whose relative an- 
gular momentum m is 1 and 3, respectively. Thus we 
can see how the ground state connects to Laughlin state 
with increasing Vi and Vs from its values of the pure 
Coulomb interaction. If there exist no phase transition, 
the ground state is in the same phase characterized by 



Laughlin state. 

In order to clarify whether there is a phase transition, 
we calculate the low energy excitations as a function of 
Vs. Since Vi of pure Coulomb interaction is already large 
in the lowest Landau level with a small probability of 
finding electron pairs whose relative angular momentum 
is 1 for z^ < 1/3, we change only V3. 

Figure 7 shows the Va-dependence of the ground state 
energy and the first excited state energy at u = 1/5 
with 12 electrons. With increasing Vs from SVs = 0, 
the ground state energy monotonically changes and the 
excitation gap continuously increases. The ground state 
correlation function at 5Vs = 0.24 shown in Fig. 6 (a) 
is circularly symmetric and it is well characterized by 
Laughlin state. Since no phase transition is detected, 
the ground state of pure Coulomb interaction belongs to 
the same phase of Laughlin state. 

When we decrease V3 from SV3 = 0, however, signifi- 
cant change occurs in the ground state. As shown in Fig. 
7, the slope of the ground state energy clearly changes 
at ^Vs = —0.04. This change in the ground state energy 
refiects the qualitative change in short range correlation 
functions. As is shown in Fig. 8, the shoulder structure 
dit X ^ 3i develops below ^Vs ~ —0.04 and it makes clear 
anisotropy around r ~ 41. The shape of the correlation 
function for SVs < —0.04 is almost the same to that re- 
alized in higher Landau levels and clearly different from 
the stripe structure of pure Coulomb interaction in the 
lowest Landau level around u = 1/5. 



3.5 Transition to Wigner crystal 

In the limit of low filling, the magnetic length i be- 
comes negligible compared with the mean distance be- 
tween the electrons. Since i is the length scale of one 
particle wave function, the electrons can be treated as 
classical point charges, and the ground state is expected 
to be Wigner crystal. In order to confirm the existence 
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Fig. 7. V3 dependence of the ground state and first excited state 
energies at u = 1/5 in the lowest Landau level. Ne = 12. 6V3 = 
corresponds to pure Coulomb interaction. 
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Fig. 8. Ground state pair correlation functions on the a;-axis at 
u = 1/5 for various 6V3. Ne = 12. SV3 = corresponds to pure 
Coulomb interaction. 



of Wigner crystal at finite we next see pair correla- 
tion functions in the ground state at low fillings. The 
ground state pair correlation functions below = 3/8 in 
Figs. 5 and 9 (a) show that the stripe structure remains 
stable down to u ~ 1/6, although it seems to approach 
hexagonal structure of Wigner crystal. 

At u = 1/6 the Wigner crystal is realized in an ex- 
cited state close to the ground state. The correlation 
function of the Wigner crystal and the energy difference 
from the ground state are shown in Fig. 9 (b) and (c). 
The energy difference almost linearly decreases with de- 
creasing u except at v = 1/5, where the ground state is 
characterized by Laughlin state. Linear extrapolation on 
u shows the crossing of the two energy levels at around 
u = 1/7. Since the momentum of the two states are 
different, we expect first order transition to the Wigner 
crystal at Uc ^ 1/7. The value of pth|e-.critical filling Uc 
is consistent with previous studies,E§E2P although exact 
diagonalizations of 6-electp3ps predict second order or 
weak first order transit ion. E2P 





0.01 




Fig. 9. (a) The ground state pair correlation function at u = 
1/6. (b) The pair correlation function of Wigner crystal at u = 
1/6. (c) Energy difference between the ground state and Wigner 
crystal at low fillings. A^e = 12. The dashed line is a guide for 
eyes. 



3. 6 Phase diagram of N = Landau level 

The present DMRG calculation shows the existence 
of four different ground states: the compressible liq- 
uid around u = 1/2, the incompressible liquid at = 
n/(2n + l) and 1/5, the weak stripe state below u ~ 0.42, 
and Wigner crystal below Uc ~ 1/7. From the above re- 
sults, the ground state phase diagram is determined as 
shown in Fig. 10. Precise critical fillings and the nature 
of the transitions are not clear for the boundary between 
the liquid states and the weak stripe state, but the tran- 
sition to Wigner crystal is expected to be first order. The 
weak stripe structure found below u ~ 0.42 has clear os- 
cillations along the stripes and this is clearly different 
from the stripe state realized in higher Landau levels. 
The ground state at v = 1/5 continuously connects to 
the Laughlin state realized in the limit of large Vi and 
Vs^ while the size of the excitation gap is vary small 
compared with the gap at u = 1/3. The pair correlation 
function at = 1/5 is strongly modified from circularly 
symmetric one, because V3 of the Coulomb interaction 
is located close to the boundary between Laughlin state 
and stripe state. 



3.7 Pairing state at half filling in N = 1 Landau level 
We next investigate the ground state in the second low- 
est Landau level. At = 1/2, fractional quantum Hall 



6 



Naokazu Shibata and Daijiro Yoshioka 



A^=0 



FQHE . FQHE FQHE^H 
1/5 Stripe 1/3 2/5 3 




-5.7 



liquid 



Fig. 10. Ground state phase diagram of the lowest Landau level 
determined by the DMRG method. 



effect has been experimentally observed. To explain 
the fractional quantization, ppwpjajpairing formation of 
electrons has been proposed.llilll3'L3^ In particular, re- 
cent studies based on exact dia^iLalizations suggest Paf- 
fiann state of triplet spin pair.Ei'LiP 

In the present study we introduce Haldane's pseudo- 
potential Vi in the second lowest Landau level. We cal- 
culate the excitation energies and the pair correlation 
functions at various Vi. In Fig. 11 we show the total 
energy of the ground state and several low energy ex- 
cited states as a function of SVi. We find two charac- 
teristic 5Vi where structure of the low energy spectrum 
changes. One is located slightly below SVi = and the 
other is around SVi = 0.03 as shown in the figure by 
arrows. Below 5Vi = the ground state is almost degen- 
erate and the pair correlation function is characterized 
by clear stripes as shown in Fig. 12 (c). For Vi > 0.03, 
the ground state and low energy excited states is almost 
independent of Vi . The correlation function is circularly 
symmetric with weak oscillations as shown in Fig. 12 (a), 
and it is similar to that observed at = 1/2 in the low- 
est Landau level, where excitation gap vanishes in the 
thermodynamic limit. Since sufficiently large SVi effec- 
tively maps the system onto the lowest Landau level, we 
expect that the ground state is compressible liquid for 
6Vi > 0.03. 

Between 6Vi = 0.0 and 0.03, the ground state is 
clearly different from the stripe state. The correlation 
function clearly decays at long distance as shown in 
Fig. 12 (b) and the degeneracy in the ground state is 
lifted with a finite excitation gap of order 0.01. Com- 
pared with the compressible liquid state at SVi > 0.03, 
the correlation function is different at short distance 
as shown in Fig. 12 (d), which shows the difference in 
the correlation function between ^Vi = 0.01 and 0.06: 
g{r)svi=o.oi - g{r)svi=o.06' In this figure, we can see 
that electrons are attracted around the origin r ^ 1.5 
consistent with the pair formation. 



3.8 Stripes away from half filling in N = 1 Landau level 
We next investigate the ground state away from half 
filling. The pair correlation functions for u = 4/9, 2/5, 
4/11, and 1/3 with 16 electrons in the unit cell are pre- 
sented in Fig. 13. These results clearly show the exis- 
tence of the stripe correlations even in the second low- 
est Landau level near half filling. In order to see de- 
tailed structure of the stripes we show the correlation 
functions on the x and y axes in Fig. 14. As is shown, 
the correlation function on the y-ax.is at u = 4/9 and 
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Fig. 11. Low energy spectrum at half filling in the second lowest 
Landau level. = 16. Inset shows the results for A^e = 24. 
5Vi = corresponds to pure Coulomb interaction. 
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Fig. 12. The ground state pair correlation functions in guiding 
center coordinates at half filling in the second lowest Landau 
level at 5Vi = 0.06 (a), 0.01 (b), and -0.02 (c). Figure (d) shows 
the difference in the correlation function between 5Vi =0.01 and 
5Vi = 0.06, g{r)sv^^Qm - g{r)sv^^o.06- Ne = 24. 



2/5 shows monotonic increase from the origin and it be- 
comes almost constant at a long distance. This is the 
same structure of the stripes observed in higher Landau 
levels near half filling as shown in Fig. 15. Compared 
with the stripes in = 2 Landau level, the amplitude of 
the oscillations on the x-axis is about 50% smaller and 
their period is 30% shorter. This means the strips in 
= 1 Landau level is likely to be disturbed by residual 
potential fiuctuations and temperatures. We think this 
may be a reason why anisotropic resistivity is not clearly 
observed in the second lowest Landau level. 

With further decreasing character of the stripes 
changes at u ~ 4/11. This is clearly shown in the corre- 
lation functions along the stripes in Fig. 14. The correla- 
tion function for < 4/11 shows weak oscillations along 
the stripes, which is similar to the stripe structure in 
the lowest Landau level. The oscillations on the stripes 
are enhanced with decreasing u and the ground state be- 
comes almost the same to that in the lowest Landau level 
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Fig. 13. The ground state pair correlation functions in guiding 
center coordinates in the second lowest Landau level near half 
fining. A^e = 16. 
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Fig. 15. The ground state pair correlation functions on the y-axis 
in guiding center coordinates in = 1, 2, and 3 Landau levels 
at z^iv = 2/5. Inset shows the results on the x-axis. A^e = 16. 
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Fig. 14. The ground state pair correlation functions on the y-axis 
in guiding center coordinates in the second lowest Landau level. 
Inset shows the results on the x-axis. Ne = 16. 



except at = 1/3. At low fillings v < 1/4, the corre- 
lation functions in 7V=0 and 1 Landau levels are almost 
identical with each other. The similarity of the ground 
state between the two Landau leiifjLs^at u = 1/4 has been 
shown by previous calculations.^ 



At u = 1/3, the ground state correlation function is 
clearly different from that of the lowest Landau level. 
Compared with Laughlin state in the lowest Landau 
level, short range correlation below r ~ 2^ is enhanced 
and clear oscillations exist at long distance as shown in 
Fig. 16. This is due to the large reduction of short range 
interaction Vi, which stabilizes Laughlin state. 

In order to see the effect of Vi on the ground state, 
we show in Fig. 17 the energies of the ground state and 
excited state as a function of Vi. The excitation gap 



Fig. 16. The ground state pair correlation functions at = 1/3 
and 1/5 in the lowest and the second lowest Landau levels. 



monotonically increases with increasing Vi from SVi = 
and no phase transition is detected. Since the ground 
state at large Vi is Laughlin state as shown in Fig. 18 (c), 
where the correlation functions at SVi = 0.08 is almost 
identical to that in the lowest Landau level, the present 
result shows the ground state at SVi = continuously 
connects to Laughlin state. However, the excitation gap 
is very small at 6Vi = compared with the gap in the 
lowest Landau level and slight decrease in Vi leads to 
the transition to stripe state as shown in Fig. 17. The 
correlation function of the stripe state at SVi = —0.02 
is shown in Fig. 18 (b), which is similar to the stripe 
state near half filling. With further decreasing Vi, first 
order transition to two-electron bubble phase occurs at 
5Vi ~ —0.04. This is consistent with the results in higher 
Landau levels, where buJ^blejState is realized at vn = 1/3 
in TV = 2 Landau leveLBMEaM 

In contrast to the ground state at vn = 1/3, the cor- 
relation function at at = 1/5 is almost the same to that 
in the lowest Landau level as shown in Fig. 16. This 
is not surprising because Vs of the second lowest Lan- 
dau level is larger than that of the lowest Landau level. 
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Fig. 19. Energy difference between Wigner crystal and the 
ground state. A^e = 12. The dashed line is a guide for eyes. 



Fig. 17. The ground state and first excited state energies in the 
second lowest Landau level at = 1/3 with A^e = 16. Inset 
shows the results for A^e = 12. 6Vi = corresponds to pure 
Coulomb interaction. 



(a) (b) (c) 

5Vi=-0.06 5V] = -0.02 5Vi=0.08 




Fig. 18. The ground state pair correlation functions in guiding 
center coordinates in the second lowest Landau level at uj^ = 1/3 
with A^e = 16. 



Since Laughlin state at j^at = 1/5 is stabilized by Vi and 
V3, large V3 compared with other Vm stabilizes Laughlin 
state similarly to the lowest Landau level, although cir- 
cularly symmetric correlation function is strongly mod- 
ified by the presence of other pseudo-potentials Vm of 
Coulomb interaction. 



3.9 Transition to Wigner crystal 

At sufficiently small un, the Wigner crystal is realized 
in the ground state even in high Landau levels, because 
electrons can be treated as point particles when mean 
distance between the electrons becomes much longer 
than typical length scale of the one-particle wave func- 
tion. In order to study the transition to Wigner crystal in 
the second lowest Landau level, we calculate the energy 
of Wigner crystal and compare with the ground state en- 
ergy. The energy difference between the two states and 



(a) (b) 

Vn = 1/4 ground state Vn = 1/4 excited state 




-10 -5 5 10 -10 -5 5 10 



X X 



Fig. 20. Pair correlation functions in guiding center coordinates 
in the second lowest Landau level, (a) The ground state at 
UN = 1/4 with 16 electrons, (b) Wigner crystal at = 1/4 
with 16 electrons obtained at a different total momentum from 
the ground state, (c) The ground state at = 1/5 with 12 
electrons, (d) Wigner crystal at = 1/5 with 12 electrons. 



the correlation functions of the guiding center are shown 
in Figs. 19 and 20, respectively. The energy difference 
between the two states monotonically decreases except at 
z/AT = 1/5, where the ground state continuously connects 
to Laughlin state. The singular behavior at un = 1/5 
is due to the decrease in the ground state energy. The 
linear extrapolation of the energy difference with respect 
to UN shows first order transition to Wigner crystal at 
UN - 1/7. 



Ground State Phase Diagram of 2D Electrons in High Magnetic Field 
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Fig. 21. Ground state phase diagram of the second lowest Landau 
level determined by the DMRG method. 



3.10 Phase diagram of N = 1 Landau level 

The ground state phase diagram of the second lowest 
Landau level is shown in Fig. 21. At half filling, the 
ground state pair correlation function suggests pairing 
state. This ground state is different from both the stripe 
state in higher Landau levels and the compressible liq- 
uid state in the lowest Landau level. With decreasing 
ujsf^ however, stripe state similar to that in higher Lan- 
dau levels is realized with small amplitude of the stripes 
(stripe I), which is 50% smaller than that in = 2 Lan- 
dau level. With further decreasing clear oscillations 
along the stripes appear below u ^ 4/11 (stripe II). This 
structure is similar to the stripes in the lowest Landau 
level. At = 1/3, the ground state belongs to the same 
phase of Laughlin state in the lowest Landau level, but 
correlation function shows clear oscillations at long dis- 
tance. This is due to the proximity to the boundary be- 
tween Laughlin state and stripe state. Below ^ 1/4, 
the correlation functions are almost identical to that in 
the lowest Landau level, and thus the ground state is 
almost the same between the two Landau levels. The 
ground state at j^at = 1/5 is expected to be in the same 
phase of Laughlin state. The first order transition to 
Wigner crystal is expected to occur at ~ 1/7. 

§4. Summary and Discussion 

In the present paper we have studied ground state 
pair correlation functions and low energy excitations of 
2D electrons in the lowest and the second lowest Lan- 
dau levels by using DMRG method. We have assumed 
completely spin polarized state and neglected Landau 
level mixing. The obtained results in the lowest Lan- 
dau level confirmed the existence of compressible liquid 
at u = 1/2, incompressible liquid at v = n/(l -H 2n) 
and 1/5, and Wigner crystal below u ~ 1/7. We have 
found new stripe state between u ~ 0.42 and 1/7, whose 
correlation function shows oscillations along the stripes, 
which is different from the stripes in higher Landau lev- 
els. In the second lowest Landau level, pairing state at 
UN = 1/2 has been confirmed. The ground state is close 
to the phase boundary between the pairing state and 
stripe state, and decrease in Vi of about 2% is enough 
for the transition to the stripe state. Between un ~ 0.47 
and 0.38 we have found stripe state similar to that in 
higher Landau levels with small amplitude of stripes. We 
have also found stripe state similar to that in the lowest 
Landau level between ~ 0.38 and 1/7. At low fillings 
< 1/3, the ground state is almost the same to that in 
the lowest Landau level. This is due to the fact that V3 



of the second lowest Landau level, which is the dominant 
interaction for un < 1/3, is slightly increased from that 
in the lowest Landau level, and this increase stabilizes 
the ground state realized in the lowest Landau level. 

We have also shown the existence of incompressible liq- 
uids at UN = 1/3 and 1/5 in the second lowest Landau 
level. The size of the gap at at = 1/3 is very small com- 
pared with that in the lowest Landau level and slight 
decrease in Vi of about 0.01 leads to the transition to 
stripe state. In actual systems, the wave function has 
finite width in the direction perpendicular to the two- 
dimensional plane. Since finite width reduces Vi from 
the value of ideal two dimensional system, incompress- 
ible liquid at un = 1/3 will not be realized in systems 
of wide quantum well. Our estimate shows the critical 
width is about for rigid density distribution in the 
direction perpendicular to the two-dimensional plane in 
square potential well. 

Recent experiment on the second lowest Landau level 
shows that there exist reentrant integerj-quantum Hall 
states around un = 0.42 and un = 0.25.li3) Our results 
suggest reentrant phase around vn = 0.42 corresponds to 
the stripe phase near half filling, although the mechanism 
of the integer quantum Hall effect is still not clear. Since 
the amplitude of the stripes is small compared with that 
of higher Landau levels, it is expected that weak random 
potentials disturb stripe correlations and insulating state 
pined by residual potential fluctuations is realized before 
forming stripe state. 

Concerning to the reentrant phase around un = 0.25, 
similar behavior is expected in the lowest Landau level, 
since the ground state obtained in the present calculation 
is almost the same between the two Landau levels. This 
means the origin of insulating behavior above u = 1/5 
in the lowest Landau level is the same to that of the 
reentrant phase around un = 0.25 in the second lowest 
Landau level. Since the energy of Wigner crystal is close 
to the ground state energy, it may be possible that ran- 
dom potentials stabilize pined Wigner crystal, which is 
the same mechanism of integer quantum Hall effect at 
small UN < 1/7. 
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